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This paper deals with the dissipative dynamics of a quantum harmonic oscillator interacting with 
a bosonic reservoir. The Master Equations based on the Rotating Wave and on the Feynman- Vernon 
system-reservoir couplings are compared highlighting differences and analogies. We discuss quanti- 
tatively and qualitatively the conditions under which the counter rotating terms can be neglected. 
By comparing the analytic solution of the heating function relative to the two different coupling 
models we conclude that, even in the weak coupling limit, the counter rotating terms give rise to 
a significant contribution in the non-Markovian short time regime. The main result of this paper 
is that such a contribution is actually experimentally measurable and thus relevant for a correct 
description of the system dynamics. 
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I. INTRODUCTION 



During the last few decades a huge deal of attention has been devoted to the study of the quantum dynamics of 
dissipative systems. The theory of open quantum systems, indeed, is essential for the understanding of a variety of 
physical phenomena in different fields of physics, such as, for example, quantum optics and solid state physics [0. 
Moreover, very recently, there has been an increasing interest in the effects of decoherence, due to the unavoidable 
coupling with external environment^, |3[ Q], on the dynamics of quasi-closed systems used for quantum computing 
and quantum information processing. The usual approach for studying decoherence and dissipation effects starts by 
prescribing a total Hamiltonian for the closed total system (system+reservoir). Then, after tracing over the reservoir 
variables and performing, if necessary, appropriate approximations, one finally derives a Master Equation ruling the 
dynamics of the dissipative quantum system. One of the most commonly done assumptions for describing open 
quantum systems is the so-called Born-Markov approximation which basically consists in neglecting memory effects of 
the reservoir, fn other words one assumes that the correlation time of the reservoir, characterizing the time scale on 
which the reservoir memory would feed back to the system, is much shorter than the typical system time scale. When 
such condition is satisfied it is possible to derive a Master Equation describing the time evolution of the dissipative 
system for times longer than the correlation time of the reservoir. Under this approximation the resulting Master 
Equation is called Markovian Master Equations and of course does not describe appropriately systems interacting 
with natural or engineered structured reservoirs, such as atoms decaying in photonic band gap materials or atom 
lasers. It has been very recently demonstrated by Ahn et al. || that non-Markovian reservoirs may be of potential 
interest for quantum information processing since a quantum system is decohered slower in a non-Markovian reservoir 
than in a Markovian one. 

In this paper we firstly derive the time-convolutionless Master Equation describing a quantum harmonic oscillator of 
frequency ujq interacting with a bosonic reservoir represented as an infinite chain of harmonic oscillators of frequencies 
{uji} |8|, ^ |l0j. The method used, exploiting a superoperatorial formalism, leads to a non-Markovian Master Equation 
spoiled of reservoir memory kernels^, In words one says that such a Master equation is local in time. Our aim is 
to analyze differences and analogies in the dynamical behaviour of this specific open system in correspondence with 
two different prefixed system-reservoir couplings. The first choise is the following: 

oo j 

Hf™ A = (sM + h.c)j = a (aft + a f i?) , (1) 

i=0 " 

usually referred to as Rotating Wave (RW) coupling. In Eq. (pi), a and bi are the annihilation operators of the system 
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and reservoir harmonic oscillators respectively and a is the adimensional coupling constant. Note that, for the sake 
of simplicity, in the paper we use adimensional position and momentum operators for the system oscillator. 

The second form of the system-reservoir interaction Hamitonian examined in this paper is the so-called Feynman- 
Vernon (FV) coupling |pTf : 

oc 

H sr = aXj2 V^I (dfbi + 9$) =aX(R + , (2) 

i=0 

where the operator X is related to the creation and annihilation operators of the quantum harmonic oscillator simply 

as 

X = i=(a + at). (3) 

While the first interaction Hamiltonian is very often used in describing quantum optics systems |Q and atom lasers 
fl3|| , the second one leads to the Master Equation for Quantum Brownian Motion fuj. Using the Hamiltonian given 
by Eq. (|l|) instead of the more general one given by Eq. (||) is usually motivated saying that the counter rotating 
terms abj and a^S^ , appearing in Eq. (^), do not conserve the total unperturbed energy and thus give a negligible 
contribution to the system dynamics in the weak coupling limit |T^| . 

The main result of this paper is that, in the non-Markovian regime, the contribution given by the counter-rotating 
terms is not negligible and experimentally measurable, also when the weak coupling limit is invoked. 

The paper is structured as follows. In Sec. II we introduce the superoperator formalism for the derivation of 
non-Markovian generalized Master Equations. In Sec. Ill we specialize the generalized Master Equations to the cases 
of Rotating Wave and Feynman- Vernon couplings and we compare them in Sec. IV. Finally in Sec. V we present 
conclusions. 

II. DERIVATION OF THE MASTER EQUATION: AN OPERATORIAL APPROACH 

Let us consider an open quantum system interacting with an environment whose physical nature needs not to be 
specified at this moment. We indicate the total Hamiltonian as follows 

H — Ho + He + aHint , (4) 

where Hq, He and Hint stand for the system, environment and interaction Hamiltonians respectively and a is the 
coupling constant. The Van Neumann-Liouville equation for the total system, in the interaction picture, is the 
following 



dp(t) a 
dt ih 



Hj(t),p(t) = ^Hf (t)p(i). (5) 



In Eq. (^), p and Hj(t) are the density matrix and the interaction Hamiltonian of the total system respectively, in 
the interaction picture, and the superoperator Hf (t) is defined as Hf (t) = [Hi(t), ■ ]. In the rest of the paper, given 
a certain operator A, we will use the following notation for the "commutator "and "anticommutator"superoperators: 

A s = [A, ■ ] A s - {A, ■ }. (6) 

In deriving the generalized Master Equation we assume that at t — -- system and environment are uncorrelated, 
that is p(0) = p(0) (^) pe(0), with p and /5g density matrices of system and environment respectively and that the 
environment is stationary, i.e. [He,Pe] = 0. 
A formal solution of Eq. (||) can be written as 

p(t) = T(t)p(0), (7) 
where the superoperator T(t) is defined as the solution of the equation: 

t(t) = ^Hf (t)T(t), (8) 
in 
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with T(0) = 1. Remembering that p(t) = 7>e {/?(£)} and that p(0) — p(0) ® Pe(0), after tracing over the environ- 
mental variables, Eq. (mj becomes 



p(t) = Tr E {T(t)p E (0)}p(Q) 

= (T(*M0)=(1+M(t))p(0), 



(9) 



where we have indicated with (T(t)) the superoperator (1 + M(i)) = 7>e {T(t)ps(0)}, acting on the space H s 
with TC S Hilbert space of the system. Differentiating now Eq. (m yields 



dp(t) 
dt 



M(*)p(0). 



Inserting in Eq. (|To|) the expression for /5(0) obtained inverting Eq. 

dp{t) 



gives 



dt 



(M(i)) [l + M(t)]-^(t) = K(i)^). 



(10) 



(11) 



In the previous equation we have defined a new superoperator K(i) acting on the space TL a ^TLg too. At this point 
it is worth spending few words on the existence of K(i), that is of the inverse superoperator [1 + M(i)] . To this 
aim, we recast K(i) in the form 



K(i) = (M(t)) [1 + M(t)]- 1 = (M(f)) i- M (t)) n 



(12) 



As discussed in [ fl5| for weak couplings such series converges at any time t. For generic coupling, however, the 
convergence radius of the series depends both on a and on t. For this reason one has always to pay special attention 
to such a problem when working for intermediate or even strong coupling regimes. 
Now, it is easy to convince oneself that a formal solution of Eq. M) may be written as 



T(t) = exp c 



^ / Hf(ti)dti 



OG 

^ (ih 



o 



t„-i 



Hf(t)Hf(ii)---Hf(t„) dh---dt n , 



(13) 



where the subscript c in the exponential stands for the Dyson chronological order, i.e. t n > i„_i • •• > t\ > t. Inserting 
such an expression into Eq. (^) with the help of Eq. (^) and collecting all the terms proportional to the same power 
in a, it is possible to demonstrate that the following expansion holds: 



(14) 



with 



k„(i) = 



tn-l 



«Hf (t)Hf (ti) • • • Hf(t„))) . c . dtx ■ ■ ■ dir. 



In the previous equation we have indicated with ((• • -)) .c. the temporal ordered cumulants [| 
report the expression of the first and second cumulants, respectively: 

«H?(i)» . c . = (Hf(t)) 
«Hf(t)Hf(^)))o. c . = (HfWHf^-^fWKHffo)). 



(15) 

As an example, we 
(16) 



The form of Eq. ( 15 ) resembles the result obtained by Van Kampen in the context of stochastic differential equations 

EE!-. . 

The origin of the expansion given by Eqs. Q14J) and (|5|) can be understood as follows. Let us write the superoperator 
K(t) defined by Eq. (13) in following symbolic form 



K(i) 



ft 



(ex Pc (^a^ Hf(ii)dti 



where 



-F[A(t)] S (iA(t))F[A(*)] 



(17) 



(18) 
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In Eq. ( [l7| ) the symbol (• • •) = tr^l ■ ■ ■ Pe} describes an operation of average over the environmental degrees of 
freedom. The expression (exp c (^a f£ Hf(£i) dt\)j} can thus be seen as the generalization, in the superoperator for- 
malism, of the concept of characteristic functional frL9| , As a consequence the superoperator In (exp c (a J* Hj(t)j ) 

is the generalization of the generator of cumulants introduced in standard textbooks. This circumstance makes it 
clear why the integrand in Eq. is called temporal ordered cumulant. In view of Eq. ([l7]) the existence problem 
of the superoperator K(t) can be traced back to the convergence of the series of cumulants in Eq.(14). 

In order to derive the explicit form of the generalized Master Equation, let us assume a bilinear interaction Hamil- 
tonian of the form: 



ffj(t) =a^2Ai{t)Ei{t) =aA(t) ■ E(t), (19) 



i=l 



where A(t) = (^Ai (i), A^if), . . . ,A n (t)j and E(t) = ^£i(£), £^(i), . . . ,E n (tfj are system and environment operators 

respectively. In the weak coupling limit we may stop the cumulant expansion given in Eq. ([l4]) to the second order in 
the coupling constant. In view of Eqs.([ll]) and (|lj) ( |l5| ) , this leads to the following Master Equation 

— = ^(tMt) 



h 2 



f [(Hf (t)Hf (h)) - (Hf (t))<Hf (h))] dh pit). (20) 
Jo 



Assuming for simplicity that the form of the environmental density matrix satisfy the condition (E) = (as for 
example in the case of a thermal reservoir), one can show that the first term of Eq.(|2^) vanishes at every time t. 
The explicit manipulation of the second term is presented in Appendix A and leads to the following final form of the 
non-Markovian generalized Master Equation: 



dp(t) 



,„ = - E (/ M"r)Af (*)Af (t - t) - mA^Af (t)A?(t - r)] dr) p(t) 

p(t)-*G(t)]p(t)=L(t)#t), (21) 



where definitions of Eq.(^) have been used. In Eq.(|T]) we have introduced the environment correlation Kij(r) 
and susceptibility /Zij(r) matrices, with r = t — t\. Such quantities, characterizing the temporal behavior of the 
environment, are defined as follows 



a 2 



KiAr)=^{{Ei(T),E j (0)}), (22) 



^Ar) = ^([Ei(r),E^0)}). (23) 

The form of Eq. (^l|) has a clear physical meaning. One can show that the superoperator D(t) is strictly connected 
with diffusion (decoherence) processes only [2(j]. The superoperator G(t), describing the dissipation processes and 
frequency renormalization, on the other hand, arises from a quantum mechanical treatment of the environment and, 
indeed, vanishes when a semiclassical description of the environment is used (see also Eq.(p3|)) |20| . 

In the next section we further carry on the calculations in order to obtain and compare the two non-Markovian 
Master Equations corresponding to the Rotating Wave and Feynman- Vernon couplings respectively. 

III. ROTATING WAVE AND FEYNMAN- VERNON COUPLINGS IN THE NON-MARKOVIAN 

REGIME 

Let us consider a quantum harmonic oscillator whose Hamiltonian is given by: 

Ho = ^ (P 2 + I 2 ) = (a f a + 1/2) , (24) 
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with wo frequency of the harmonic oscillator. The system interacts with a bosonic reservoir at T temperature of 
Hamiltonian 



H E = hJ2^i (SJS* + 1/2) 



(25) 



i=0 



with u>{ frequencies of the reservoir oscillators. 



A. Feynman- Vernon coupling 



Let us begin discussing the Feynman- Vernon interaction Hamiltonian, given by Eq.pl). Such a coupling is of the 



form of Eq.(19) where, in the interaction picture, A(t) = X(t) and E(t) = R(t) + iv(t). Our aim is to manipulate 
Eq.(pT|) in order to obtain the specific non-Markovian generalized Master Equation appropriate for our system. In 
this section we will sketch the main steps of the derivation. More details can be found in Appendix B. 

First of all let us write the Master Equation given in Eq.(^l|) in the Schrodinger picture. Introducing the superop- 
erator 



T (t) = exp 



ih 



-mt 



(26) 



with Ho given by Eq.(|2^), and transforming in the Schrodinger picture the superoperator L(t) defined in Eq.(pT|) 

L s (() = T ( ( )L(i)T ^(f), (27) 

our generalized Master Equation becomes 



dpsjt) 
dt 



-n$ -B s (t)+iG s (t) 



Ps(t), 



(28) 



with ps density matrix of the harmonic oscillator in the Schrodinger picture. In Appendix B we show that the 
superoperators Ds{t) and Gg(t) can be recast in the form 



T> s (t) = [ k(t)X s (cosw tX s - sinw rP s ) dr, 
Jo 



(29a) 



G s (t) = [ p(r)X. s (cosw tX s - sinw rP s ) dr, 
Jo 



where P and P are the "commutator" and "anticommutator" superoperators associated to the operator 

*-£(«•-«), 



(29b) 



(30) 



and k(t) = «i,i(r) and /x(r) = /ii,i(r) are defined by Eqs. (E2)-(E3[). Inserting Eqs.(29a)-(29b) into Eq.(|28|) we get 



dps(t) 1 „ s 
— ±lfi 



A(i)(X s ) 2 - n(t)X s P s - -r(t)(X 2 ) s + i 7 (i)X s P s 



dt ih 

The time dependent coefficients appearing in the previous equation are defined as follows 



ps(t). 



A(t) 

7(*) 
U(t) 



k(t) cos(uJoT)dT, 



/Lt(r) sui(ujoT)dT, 



k(t) sin(wor)cir, 



r(t) = 2 / p(r) cos(w r)(iT. 
'0 



(31) 

(32) 
(33) 
(34) 
(35) 
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From the form of Eq. ( |3l| ) , and remembering that Hg may be written as 

H£ = i[(P 2 )* + (X 2 f], 



(36) 



it is not difficult to convince oneself that the term having coefficient r(t) gives a renormalization of the oscillator 
frequency. 

As usually done in standard textbooks Jl2| , this term can be included in the definition of u>q. In the following we 
neglect such term since it is possible to prove that such an approximation is always justified in the weak coupling 
regime a -C 1, provided that the reservoir frequency cut-off remains finite. 

Under these conditions, the Master Equation, in the interaction picture with respect to H , assumes the form 



dp® 
dt 



[A(t)(X s ) 2 (t) - U(t)X s (t)P s (t) + z 7 (i)X*(i)P s (t)] p(t). 



The time dependent superoperators appearing in Eq.(|37|) are those related to the operators 

X(t) = X cos(ui t) + Psm(u> Q t), 
P(t) = Pcos(w *) - Xsm(u t). 



(37) 



(38) 
(39) 



Eq. (|37| ) can be exactly solved in an operatorial way and the solution has an operatorial form |22] . This fact may be 
exploited to fully disclose both the short time non-Markovian and the asymptotic Markovian behaviors characterizing 
the dynamics of the system, as we will see in Section IV. 

B. Rotating Wave coupling 

The generalized Master Equation correspondent to the interaction Hamintonian given by Eq. (|f|) , derived following 
the same procedure presented in the previous subsection, is (see also Appendix C) 



dm 

dt 



A (t)a' a' 



a' a s a s a f 



„RWA 



(t) (at' 



„RWA( 

— v _. _ _ _ j , 2 . 

The time dependent coefficients appearing in this equation are defined as follows 

ft 



a s +a s at 2 



A RWA (i) 
7 RWA (i) 



« RWA (r)dr, 



M r A (r)dr 



r KWA (i) - 2 / M f WA (r)dr, 
Jo 



p(t). 



(40) 

(41) 
(42) 
(43) 



where 



K *WA (r) = |_ ({ ^ (r) ^t (0 ) } + {i ?t (T) ^ (0)}); 



(44a) 



(44b) 



pf WA (r) = 1 ^([R( t ),RH0)] - [&( T ),R{0)]), 



(44c) 



and R(t) is the operator (see Appendix B) 



R(t) = *£NT9ib 



-u> )t 



(45) 
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In the Schrodinger picture Eq. ( |40| ) takes the form 



df>s{t) 
dt 



-H S + zr RWA (i) (at a ) S - A RWA (t)at S a s - 1 ^ (W - a s a t S 

in ' } v 



hit) (46) 



The term proportional to r RWA (i) gives rise to a renormalization of the oscillator frequency as for the Feynman- 
Vernon case, described in the previous subsection. Therefore, proceeding with the same considerations and passing 
to the interaction picture, we obtain the following generalized Master Equation for the system 

dp A R WA(, )+7 RWA W 

= - [a'ap - 2apa> + pa'a\ 

_ A (t)_j W [aaV-2&tpa + paa t ] . (47) 

Let us note that this Master Equation, differently from the one obtained for the FV coupling (see Eq. (|3^)), is in the 
Lindblad form as far as the time dependent coefficients A RWA (t) ±7 RWA (i) are positive. This is usually the case for 
typical reservoir spectra and parameters, as we have discussed in |23|. 



IV. COMPARISON BETWEEN THE RW AND THE FEYNMAN- VERNON COUPLING MODELS 

In the previous section we have seen that starting from a FV coupling or from a RW coupling of an harmonic 
oscillator with a thermal reservoir it is possible to obtain a generalized Master Equation local in time describing the 
dynamics of the oscillator. This fact is not surprising. Indeed, as underlined by Paz and Zurek in "perturbative 
Master Equations can always be shown to be local in time ". It is worth noting that, as far as the FV interaction 
model is concerned, an exact Master Equation, valid for every value of the coupling strength, has been derived p3] . 

The Master Equations we have derived in the paper are based on the weak coupling assumption but do not rely 
on the Born-Markov approximation so we are able to examine the non-Markovian short time behavior of the system 
under study. In addition such equations of course describe the correct Markovian long time asymptotic behavior |p3| . 

The different structure of the two Master Equations given by Eq. (|37]) and Eq. (^), traceable back to the two 
different coupling Hamiltonians, are responsible for the occurrence of some physically transparent changes in the 
oscillator dynamics, more marked in short time regime. 

To better understand the physical origin of such differences let us have a closer look at the two interaction Hamil- 
tonians: 

ff RWA = a (ai? f + a f i?) (48a) 



Hl v = a 



oi? t + a^R) + (aR + a f 7? f 



(48b) 



We take advantage of a pictorial representation of the four different interaction terms appearing in Hamiltonian 
( ggg ) (see fig |). 



The events represented in the first two diagrams (a,b) are processes of absorption or emission in which energy is 
conserved. The last two diagrams (c,d), on the contrary, describe events not corresponding to real absorption and 
emission processes. For this reason such processes are called virtual processes. In the second order in perturbation 
theory both the two real and virtual processes combine to give rise to real processes hereafter called alpha and beta 
processes respectively (see fig.||). 

Thus when we use the Feynman- Vernon coupling instead of the Rotating Wave one, the channels through which the 
oscillator exchanges energy with the reservoir are doubled. The asymptotic long time behavior describes, of course, 
thermalization in both cases. 

These particular features give rise to different predictions of the short time behavior of physical quantities, such as 
for example the mean number of quanta (n){t) of the system oscillator (heating function), depending on which of the 
two system-reservoir coupling models is used. We show in the following that such different behaviors are, in principle, 
experimentally observable and thus relevant for the correct description of the complete dynamics of the system. 

Let us consider, as initial state of the system, the vacuum. It is well known that, in this case, due to the interaction 
with the thermal reservoir at T temperature, the system experiences heating processes leading to thermalization. In 



a) 

a R 



c) 

aR f 



b) 



aR 




d) 



aR 




FIG. 1: There are four distinct terms in the Hamiltonian given by Eq. (48t). The events represented in the first two diagrams (a,b) 



correspond to real processes. The last two diagrams (c.d), instead, describe events corresponding to virtual processes 





Process alpha 



Process beta 



FIG. 2: Two real processes combine to give rise to the real process alpha. Two virtual processes combine to give rise to the real process 
beta. 



reference j23| it has been shown that, for the FV coupling, the non-Markovian time evolution of (h)(t), in the weak 
coupling limit, is given by p2, G3fl 



(n)(t) 



2a 2 



uj\g(uj)\ 2 (n(ijj) + — )du> 



t z 



(49) 



where n(uj) is the mean number of reservoir excitations at T temperature and g(u>) is the reservoir spectral density. 
For the RW coupling, similar calculations yields the following expression 



\RWA 



(*) 



t«UJ~ 



Lj\g{bj)\ 2 n(oj)dijj 



l t 2 



(50) 



Comparing these last two equations one sees immediately that, for short time intervals, (n)(i) » 2(n) RWA (i), 
meaning that the system-reservoir FV coupling model predicts an initial heating of the system faster than the one 
predicted by the RW coupling model. This fact can be easily traced back to the doubling of channels for energy 
exchange illustrated in Fig. (1). 

On the other hand, in the long time asymptotic limit, (n)(t) and (n) RWA (i) have, as expected, the same temporal 
behavior 22, p3| : 



(n)(t » uo 1 ) = (n) HWA (t » lUq 1 ) ~ n(w )(l - exp[-ira 2 uJo\g(uJ )\ 2 t}), 



(51) 



since, due to the time-energy uncertainty principle, for long times t, (3 processes (see Figure ^) are very unlikely to 
happen in the weak coupling regime. 

Summing up the two system-reservoir coupling models under scrutiny predict the same asymptotic long time 
behavior for the observable (h)(t) but different non-Markovian short time behaviors. It is worth noting that, once 
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known the system and reservoir parameters, the only phenomenologic constant is the coupling constant a. Such 
quantity is usually estimated from the experiments p5| . If we now assume that experiments may be performed in all 
the relevant time scale, that is both in the asymptotic long time regime and in the non-Markovian short time regime, 
one can use the value of the coupling constant experimentally measured in the asymptotic long time regime (see Eq. 
(pl|)) to verify if the correct short time behavior is actually the one predicted by Eq. ( [f9|) (FV coupling) or the one 
given by Eq. ( |50| ) (RW coupling). In fact, one would expect that, since the complete Feynman- Vernon coupling is 
more general than the RW coupling, it is also more fundamental and thus it should give the correct description of the 
dynamics of the system. 



V. THE RWA IN THE FEYNMAN- VERNON MODEL: COMPARISON WITH THE RW MODEL 



Let us now consider again the final form of the generalized Master Equation, given by Eq. (37) with Eqs. 
(|39|), derived for the FV coupling. To further simplify the calculation one could think to perform a Rotating Wave 
Approximation (RWA) averaging on an interval At the rapidly oscillating tri gon ometric functions appearing in Eq. (|37~ 
through Eqs. j3q)-(|3S|) . Under such conditions, that is for 2ui At 1, Eq. ( |37[ ) assumes the form 



dp(t) 
dt 



A(t) 



((X s ) 2 + (P s ) 2 ) +i 



M-t) 



X b P^ P^X 2 



p{t). 



(52) 



Having in mind Eqs.(||) and (^), after some straightforward calculations the following Master Equation is obtained 



dp 

dt 



— — — [a/ap ~ 2apa/ + pa? a] 



A(t) - 7 (t) 



\da) p — 2a? pa + pad/] 



(53) 



53) is in the Lindblad form as far as 
221. 



with A(f) and j(t) defined by Eqs.(p2[)-(p3|). It is important to note that Eqs. 
the sum and difference coefficients (Am — j(t)) and (A(t) + 7(f)) are positive 

Another interesting feature of Eq. (|5^) is that it has the same structure of the Master Equation obtained starting 
from the RW coupling (see Eq. (47)). Indeed, one sees immediately that the difference between the Master Equation 
obtained starting form the FV coupling and performing after the RWA, and that one obtained starting from the 
RW coupling relies only on the time dependent coefficients of the ME. Let us have a closer look at the form of such 
coefficients. In the limit of continuous modes they are written as: 



A(t) = 2a 2 



o 



u\g{uj)\ 2 (n(u)) 



1 . 



cos[wr] cos[u>QT\dwdT, 



(54a) 



ui\g(<jj)\ (n(w) + -)cos[(w — u>o)T]du>dT, 
o 2 



(54b) 
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7(t) = 2a 2 / / —\g(uj)\ 2 sin[LUT}siii[ujoT}dujdT, (55a) 



7 KWA (t)=a 2 / / -\g{io)\ 2 cos[{uj - uj Q )T]dwdT, (55b) 



where n(w) = (earp[fe] — l) is the number of reservoir excitations at T temperature. In the following we assume 



an Ohmic environment characterized by a reservoir spectral density having frequency cut-off uj c , as for example the 
Drude spectral density 



\^ 2 = \^h • (56) 

A noticeable difference between the A(t)(j(t)) and A RWA (t)(7 RWA (<)) coefficients is that in the last one the anti- 
resonant term cos[(cj + ojq)t] is absent. Such a circumstance leads to distinguishable short time behaviors of the FV 
and RW coefficients. 

It is indeed possible to prove that in A(t), for t <C alpha and beta processes give rise to the same contributions 
linear in t so that A(i) » 2A RWA (i). 

As far as j(t) is concerned, on the contrary, the same processes cancel each other at the first order in t in such a 
way that, for t <C w^ 1 , j(t) oc t 3 whereas 7 RWA (i) oc t. 

In the asymptotic Markovian long time regime we have, as expected, that A(t 3> ui^ 1 ) — A RWA (t ^> w" 1 ) and 
j(t > cj" 1 ) ~ 7 RWA (i > cj- 1 ), as shown in Fig. | 

At this point it is worth making some considerations on the validity of the RWA performed to derive Eq.(|53|). As we 
have already said at the beginning of this section, the RWA consists in neglecting terms oscillating at the frequency 
2ujq. In other words performing the RWA amounts at looking at the course-grained structure of the dynamics of the 
systems. For this reason we cannot describe correctly the dynamical features in a time interval such as At < loq 1 . 
Very often one deals with situations in which the characteristic frequency of the system loq is smaller or much smaller 
than the reservoir frequency cut lu c . Under this circumstances, normally, we cannot rely on the short time expressions 
of the FV coefficients A(t) and j(t) since they are valid for times t -C w^ 1 <C ^q 1 - However, there are two cases 
in which one can use the Master Equation given by Eq. (|53|) to describe correctly the non-Markovian short time 
behavior of the system: 

1. whenever one wants to look at situations in which loq > lu c , as discussed for example in Jl^, [2^] ; 

2. whenever we are interested in the mean value of a certain class of observables, like for instance the number 
operator (ata)(t) (see |p2|). 



In this last case, indeed, it has been shown |22 that, in the weak coupling limit, it is equivalent to use the solution of 
the the Master Equations ( |37| ) or (||^) since they lead to the the same analytic expressions for the expectation value 
of the observable of the class before mentioned. 

VI. SUMMARY AND CONCLUSIONS 

We have described a procedure, based on superoperator formalism, to derive, in the weak coupling limit, non- 
Markovian generalized Master Equations local in time. Such a method is equivalent to the time-convolutionless 
projection operator technique in the sense that it leads to the same generalized Master Equation. Wc apply this 
procedure to derive the Master Equation for a specific system, namely a quantum harmonic oscillator coupled to a 
thermal reservoir at T temperature. We compare two different microscopic system-reservoir coupling models: the 
Feynman- Vernon and the rotating wave couplings. Both couplings are bilinear, but the first one is more general 
and thus, in this sense, it is more fundamental. Very often however, in quantum optics systems, the Rotating Wave 
coupling is used because the counter rotating terms not conserving the unperturbed energy cannot contribute to the 
system dynamics . The main result of our paper is to establish under which conditions such a claim is effectively 
correct. By comparing the analytic solutions of the heating function relative to the two different coupling models 
(FV and RW couplings) we conclude that, even in the weak coupling limit, the counter rotating terms give indeed 
a significant contribution in the non-Markovian short time regime. Such a contribution is actually experimentally 
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measurable, provided that one can perform experiments in all the time scale relevant for the system dynamics. To this 
purpose it is worth noting that in the context of trapped ions experiments have been performed in which the system 
(single harmonic oscillator) is first cooled down to its zero point energy and then coupled to a properly engineered 
reservoir |p5[ . We note that, in such experiments, it is possible not only to choose at will the reservoir parameters, 
but also to engineer the coupling and control the coupling strength. Therefore, the great experimental advances of 
the trapped ion techniques could make it possible to perform an experiment aimed at proving the relevant role, in 
the short time dynamics, of the usually neglected counter rotating terms. 

One of the reasons for which one usually prefers to work with master equations derived starting from the RW 
coupling model is related to the fact that the resulting Master Equation, in the Born-Markov approximation, is in the 
Lindblad form differently from the case in which the Feynman- Vernon coupling is assumed (see Master Equation for 
Brownian motion). We have demonstrated here that also the non-Mar kovian Master Equation obtained starting from 
the RW coupling is in the Lindblad form, for some value of the relevant system and reservoir parameter. Moreover, 
by looking at the analytic expression of the time dependent coefficients of our non-Markovian generalized Master 
Equations one can infer the conditions under which one passes from Lindblad to non Lindblad Master Equations. 
Remembering that the Master Equations given by Eqs. ([l7|) and ([53| ) are of Lindblad type when their time dependent 
coefficients are positive, indeed, it is not difficult to convince oneself that such conditions are simply related to the 
change of the sign of the coefficients. Therefore the form of the RW Master Equations derived in this paper allows 
to study the border separating two very different physical regimes characterized by very different system dynamics 
fl5f and, for this reason, makes it possible to gain more insight in the fundamental dissipative processes of one of the 
most extensively studied physical systems: a harmonic oscillator coupled to a thermal reservoir. 

Another new result we have obtained in this paper stems from the comparison between the master equations derived 
in the following two cases: 

1) Feynman- Vernon system reservoir coupling followed by the RWA performed after tracing over the reservoir 
degrees of freedom; 

2) Rotating Wave system reservoir coupling. 

Stated another way we look at the differences in the system dynamics arising from the two following approximations 
respectively: 

1) average over rapidly oscillating terms after tracing over the reservoir variables; 

2) neglecting the counter rotating terms in the initial microscopic coupling model. 

We have shown that the Master Equation obtained from the Feynman- Vernon coupling, after performing the RWA, 
is of the Lindblad type and it actually has the same structure of the RW Master Equation, with different time 
dependent coefficients. We have demonstrated that these two different approximations lead to different short time 
behaviors, while in the asymptotic long time Markovian regime the two correspondent Master Equations do coincide. 
However we have proved that performing the RWA after tracing over the reservoir variables is a less restrictive 
approximation than starting with the RW coupling model. Indeed, differently from the RW Master Equation, the 
Feynman- Vernon one + RWA, takes into account the virtual photon exchanges relevant in the short time dynamics 
and thus it predicts the correct non-Markovian short time behavior, provided that t 3> 0Jq . 



One of the authors (S.M.) acknowledges financial support from Finanziamento Progetto Giovani Ricercatori anno 
1999, Comitato 02. 



In this Appendix we sketch the derivation of Eq. ( |2l| ) from Eq. ( |20| ) . Let us consider an interaction Hamiltonian 
with the form 
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APPENDIX A: THE GENERAL MASTER EQUATION 



Hi{t) = aA t (t)E t (t) = al(t) ■ E(t), 



(Al) 



where for simplicity we have used Einstein notation. Using some algebraic properties of the superoperators one can 
show that if [Ai, Ei] = then 




(A2) 
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Exploiting the properties of the trace and the assumption (E) = one gets 

tr E {E?p E (0)} = 0, tr E {Efp E (0)} = 0. 

Consequently 

(Hj S (t)) = tr E {Efp E (0)} = \ ((Ef (t))Af (t) + (Ef (f))A?(f)) = 0. 
In the same manner it is not difficult to show that 

(Hi S WH, s (^)) = \ «Ef (t)E= (h))Af (t)Af (h) + (Ef(t)Ef (t x ))Af (t)A?(ti)) 
where the equalities 

(Ef(t)E 3 s (t 1 )) = (Ef(t)Ef(i 1 )>^0, 
have been used. After some algebraic manipulation one gets 

(Ef mf( tl )) = 2{{E?(t - h),Ef(0)}), 



(Ef(t)Ef(t 1 ))=2([4 s (t-t 1 )^f(0)]), 



(A3) 
(A4) 

(A5) 

(A6) 

(A7) 
(A8) 



where the square an curl brackets indicate the commutator and anti-commutator respectively. 

Substituting these expressions into Eq.(p0|) and using the definitions of the correlation and susceptibility matrices 
one obtains the general Master Equation given by (^l|). 



APPENDIX B: DERIVATION OF FV MASTER EQUATION 

In this appendix we present the superoperatorial mathematical properties allowing to derive the final form of the FV 
Master Equation, given by Eq. (37) discussed in this paper. First of all let us consider the following superoperatorial 
relations 



[A S ,B S ] = [A S ,B S ] = [A,B] S 
[A S ,B S ] = [A,Bf. 

Such these equations and having in mind the Baker-Hausdorff formula one gets 



A s(s) (t) = exp [iB s t] A s(s) exp [-iB s t] = (exp 



iB s t 



A exp 



-iB s t 



S(S) 



(Mtj) 



S(S) 



(Bla) 
(Bib) 



(B2) 



The previous relation says that the time evolution of superoperators is ruled by equations analogue to those of the 
operators in Dirac's formalism. Then, specifying Eq. ( ^l| ) to the system under scrutiny and using Eqs. ( p^ ) and ( pgj ) 
yields 

D s (t) = f K (T)T Q (t)X s (t)X. s (t-T)T Q 1 (t)dr= f K (r)X s X s (-r)dr, 
Jo Jo 

= [ k(t)X. s (coscj tX s - sincjoTP s ) dr 
Jo 



(B3a) 



G s (t) = /V(T)ToWX s (t)X s (t-r)T- 1 Wdr: 
Jo 

= / m( t ) xS (cos cj tX s - sin w rP S ) dr. 
Jo 



/x(r)X 6 X 2j (-r)(ir 



(B3b) 
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APPENDIX C: DERIVATION OF THE RW MASTER EQUATION 

In this Appendix we underline the essential steps in the derivation of the RW Master Equation given by Eq. (|4 
Let us note that, in interaction picture, the Hamiltonian given by Eq. (nj) assumes the form 



oo j — 
a — n * 



+ h.c). (CI) 



From a mathematical point of view it is convenient to associate all the time dependent phase factors to the reservoir 
operators as follows 



££ WA (i) = a (aBl(t) + a*R(t)) , (C2) 

with 

OO I — 

^)=E fi \6«^" (Ui "" o)i ( C3 ) 

For the system here considered, thus, the operators appearing in the bilinear form defined in Eq. ( JTg| ) are A(t) = 

A= (a, at) and E{t) = (t), R(t)). 

Exploiting the properties of superoperators given by Eq. (Bl), with some algebraic manipulation, the Master 
Equation given by Eq.(^lj) can be recast in the form given by Eq. (|4^). Finally we write such a Master Equation in 
the Schrodinger picture exploiting of the following property 

T (t) atV( s > To Ht) = (at e — °*) s (a e^" 1 ) 3 ^ = a* V< E >, (C4) 



with To(t) defined by Eq. (|2q). Concluding, we note that, the superoperator proportional to r RWA (i), appearing in 
Eq.(pT|), can be recast in the form 



a 6 at^) = (fit a)". (C5) 



Thus the corresponding term in the RW Master Equation is a frequency renormalization term. Neglecting this term 
and going back to the interaction picture one gets the final form of the Master Equation, given by Eq. (f47|). 
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